Effective Breather Trapping Mechanism for DNA Transcription 
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Collective coordinate and direct numerical integration methods are applied to the analysis of 
a one- dimensional DNA model. A modification of the coupling constant in an extended region is 
found to be less selective towards the breather it can trap than an isolated impurity. Therefore it 
provides a possible physical mechanism for the effect of an enzyme on DNA transcription. 
PACS: 87.10.+e, 03.40.Kf 

I. INTRODUCTION 



The first step of the transcription of deoxyribonucleic acid (DNA) is a local opening of the double helix which 
extends over about 20 base pairs. Such local unwindings of the helix can be obtained by heating DNA to about 
70° C. But in the life of an organism they must occur at physiological temperature. This is achieved by the action 
of an enzyme Q. However one may wonder how this can be possible since, whatever its origin, the local opening 
| requires the breaking of the same number of hydrogen bonds, hence the same amount of energy, and the enzyme does 
. not bring in energy. However, under normal physiological conditions there are thermal fluctuations along the DNA 
f — ' chain. They can be weakly localized by nonlinear effects to generate what biologists call the "breathing of DNA" . But 
their intensity is not high enough to open the double helix over many base pairs. A possible pathway to the opening 
would be to collect the thermal energy that is present along the molecule. This could be the role of the enzyme. From 
a physicist point of view, the effect of an enzyme can be considered as a perturbation to the DNA lattice. 

Recently Forinash et al considered the interaction between a mass impurity on a DNA chain, and thermal nonlinear 
waves described as breathers traveling along the chain. [Q They found that the impurity is selective toward the 
. breather it can trap. Although this is a first indication that a defect can contribute to localize energy in a nonlinear 
chain, it does not appear to be a good model for the action of an enzyme because, with such a localized defect, only 
some predefined frequencies of the thermal fluctuations would contribute to bring in the energy. Therefore one may 
I ■ ask whether there exist any other mechanism more efficient to trap energy. 

One learns from biological studies that some proteins, make contact with DNA at multiple sites J^,Q]. Moreover the 
transcription enzyme actually bends DNA toward itself, ft has the effect not only to modify the mass at some sites 
but also to modify the coupling constants along the strands. The bases which are inside the bend are brought closer to 
each other while the ones which are outside are moved farther apart. Although the variation of the distances between 
neighboring bases may be rather small, it can have a large effect because the interaction between bases is due to the 
overlap of tt electrons over the whole surface of the planar bases. We examine in this paper whether the interaction 
of the enzyme with more than one site might be more efficient for trapping breathers than isolated impurities by 
studying the effect of an extended modification of the coupling along the DNA chain. 

The effect of bending and twisting to modify the elasticity of DNA has been considered previously by Barkley and 
Zimm ||, and by Marko and Siggia || but they did not study the consequences of base pair opening. Salerno 
considered the dynamical properties of a DNA promoter which has some similarities with our problem because we 
treat here the enzyme as an inhomogeneity due to an external effect while he considers inhomogeneities from the DNA 
composition itself. However he was interested by kinks while we study breathing modes. In a more abstract level we 
are investigating here a nonlinear model, with an "extended defect" , and we try to understand the interplay between 
nonlinearity and disorder. In the harmonic case, a one dimensional chain with isolated defects has been considered 
before by Montroll and Potts ||. However, besides the introduction of nonlinearity, one should also notice that for 
the type of extended defect that we consider, there is no evanescent local mode which would couple to a breather as 
in the case considered by Forinash et al, so that the mechanism for energy localization must be different. 
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II. DNA LATTICE MODEL 



If one neglects the small longitudinal motion and concentrates on the stretching of the base pairs, DNA can be 
described by a simple one-dimensional model |)| which consists of an array of harmonically coupled particles subjected 
to a Morse potential. Such a model is sufficient to provide a good qualitative description of the thermal denaturation 
of the molecule . If one treats a chain with inhomogeneous coupling, the equations of motion read 

d 2 Y 

m -dfT - K n+ i{Y n+1 - Y n ) + K n (Y n - y„_i) 

-2Dae- aY -(e' aY " - 1) = , (1) 

in which a and D are parameters for the Morse potential jyj , which have dimensions of inverse length and energy 
respectively, and n is the site index. Fig. |l| shows the geometry and the coordinate used. It is convenient for the 
analytical calculations to transform these equations into a dimensionless form by defining the following dimensionless 
variables: 

Vn = aY n , (2) 



t - J^T, (3) 
V m 



The equations become 

d 2 y 



g2 k n+1 {y n+1 - y n ) + k n (y„ - y n -i) 

-2e- yn {e- yn -1) = 0. (5) 

One notices that the last set of equations contain only one parameter, the coupling constant. In order to represent 
the perturbation due to the enzyme, one could imagine to modify locally any of the parameters of Eq. ([!]), but it is 
likely that the presence of an enzyme will affect the coupling constant through the bending of the molecule. Moreover 
previous studies of the role of disorder on the dynamics of the DNA model |l2j have shown that the formation of 
open regions in the model are much more sensitive to modulations of the coupling constant than to changes in other 
parameters. Therefore we only consider here an extended perturbation of the coupling constant. An additional 
possibility to model the enzyme specificity is however examined in the discussion. 

Since we do not know how to solve in the discrete case, we transform the set of equations, Eqs. (|^), into the 
corresponding continuous PDE. In the continuum limit, with a Taylor expansion in the potential term which assumes 
small amplitude oscillation, Eq. (|5|) becomes, 

in which d is the lattice spacing and k\ is a space dependent coupling constant. We set d equal to unity in the 
following calculations. 



III. THE PERTURBED NONLINEAR SCHRODINGER EQUATION 

Equation (|^) can be transformed into a perturbed Nonlinear Schrodinger equation by a multiple-scale expansion 
L3|,0. Assuming that the amplitude of the thermal oscillation is small, y w e<p, we perform the expansion 

5(e 3 ) , (7) 
+ 0(e 3 ), (8) 

-e 2 + 0(e 3 ). (9) 

> 

Moreover we assume a modulation of the coupling constant of the order of e, i.e. 
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— « — e (10) 
<9a; 9a;i 



Equating like powers of e yields a sequence of equations, in ascending powers of e: 

(^i + o 92j? o ) _ 9fci_^o _ fco( ^i + 2 d2p o ) 
di 2 , dt dti dx\ dx Q dx 2 , dx dxi 

+2(F 1 - ^F 2 ) = , (12) 
and higher order equations , 

in which ko is the unperturbed coupling constant. Solving for equations in each order of e sequentially one obtains, 

F = u[x 1 ,X2,t 1 ,t 2 )e i(qXo - ut ^ +c.c. , (13) 
3 In 2 

P. - -U 2 -l 2i(?xo-w*o) 

+c.c. , (14) 

and the dispersion relation: lo 2 = uj 2 + k$q 2 , with to 2 = 2. From the vanishing of the secular equation at q = one 
obtains the perturbed Nonlinear Schrodinger equation (NLS) at order e 2 : 

„. du dki du , d 2 u „ . , 2 

2 ^ + ^^ + fcl ^f + 8uH = - (15) 

We can further rescale the equation into a standard form: Defining the following new dimensionless variables, 

k(x) = ^ - 1 , (16) 

(17) 

^ 2 




2/c 

* = , (19) 

with k being the normalized deviation coupling in the vicinity of the enzyme, one obtains the following perturbed 
dimensionless NLS, 

1 19- 
iui + -u$x + u\u\ 2 + (ktix) = , (20) 

and the corresponding Lagrangian density, 



i , * + „ 1 ?i . - |2 , 1 1 „-|4 



A = - M|) - -(1 + fc)!^! 2 + -|zi| 4 . (21) 



2 V 1 " 2 V " 1 2 
In the following section we drop the * for nomenclature simplicity. 



IV. ONE SOLITON COLLECTIVE COORDINATE ANALYSIS. 



The collective coordinate method, which is a particle description of the soliton in contrast to the field description 
given by the Lagrangian, provides a good way to study the influence of a perturbation on a soliton. The spirit is the 
same as using the center of mass to analyze the behavior of a system of particles. 

Without the perturbating term in Eq. (^p|), one has a breather solution 
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u(x, i) = 77sech[ry(a: - u e t)]e iu «(*- u "*) + c.c. , (22) 



in which r\ — \/(u 2 — 2u e u c )/(2PQ), where u e is the envelope velocity, u c the carrier velocity, and P = 1/2, Q — 1 are 
coefficients of the second space derivative and the nonlinear terms in Eq. ( pp| ) respectively. In view of this solution, 
we use an ansatz for the collective coordinate analysis 

u(x,t)=rj sech^x - Q e^* 4 *^ , (23) 

where the parameters rj, £, <f>, £ are functions of t. For an unperturbed system this implies the following relations 
between the parameters: 

V = V^i ~ 2 u e ii c , (24) 
C = u e r)t, (25) 
Z = u e , (26) 
= —u e u c t . (27) 

At t = 0, C = and = and there are only two parameters left, which is consistent with Eq. (|2^), because the NLS 
breather is a two-parameter solution. Even when the breather is far away from the defect, because the ansatz extends 
to infinity and always feels the defect, we do not expect these relations to hold for a perturbed system. Hence in what 
follows we examine the whole four-parameter space for the equations of motion. 

Introducing this ansatz into the Lagrangian density, Eq. (glj), and integrating over space, one obtains an effective 
Lagrangian, 

L = -2r)<p t - 2C6 + 4r " £ »7 " o / k\u x \ 2 dx . (28) 



and the corresponding Hamiltonian 
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„3 1 r+oo 

H=-±-+e V +- k\u x \ 2 dx, (29) 



which contains no momentum term. 

At this point, we must specify an expression for k(x) to proceed. For algebraic convenience let us choose 

k = k[Q{x + I) - Q(x - l)\ , (30) 

in which is the Heaviside step function and I is the half-length of the defect. This form of k violates Eq. (|l0|), however 
previous works showed that the collective coordinate results are generally robust for the treatment of dynamics in the 
presence of perturbation [fi5f , therefore we can expect to get results which are at least qualitatively correct in spite of 
this rather crude approximation. Moreover we shall check them against full numerical simulations in the next section. 
Introducing the following abbreviated notation: 

T+ = tanh(? ? ; + C) , (31) 
T_ = t&nh(r)l - () , (32) 
S+ = scch(r)l + C) , (33) 
S- = scch^Z - C) , (34) 

one obtains 

r + oo 

k\u x \ 2 dx = -(Tl +T!)r, 3 + k(T + +T_)e 2 ry , (35) 

> 

which characterizes the effect of the defect and decays fast towards zero soon outside of the impurity region, and the 
equations of motion: 

* = y - y - ~ A (T + + T-)e - y (S 2 + + S 2 _)ev 

-^(S 2 + T 2 + + S 2 _T 2 _)rf- ^(Tl+T^ 2 , (36) 
6 = --AS 2 + T 2 + S 2 _T 2 W - -(S\ S 2 _)ev, (37) 



Ct=tV + %(T + +T-)tTi, (38) 
Vt = • (39) 
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As expected, far away from the defect, i.e. when S and T vanish, one recovers the usual relations for the NLS 
equation because in this case £t = so that £ is a constant that we can denote by u e . Then Q = £77 gives £ = M e ryt as 
expected, and 4> t = (?7 2 — £ 2 )/2 gives </> = —u e u c t if u c is defined by Eq. ( |24| ) for 77. 

In the presence of the defect, the set of nonlinear differential equations for the collective variables cannot be 
integrated analytically. It is however much simpler than the full set of discrete equations since it contains only 3 
equations. It can be integrated by a fourth-order Runge-Kutta method. One can however make general remarks on 
the properties of the solution before resorting to numerical calculations. The soliton described by the ansatz is an 
unbreakable entity and moreover the energy given by Eq. (|29| ) is conserved even when a potential well is encountered. 
Therefore when the soliton reaches a defect, it may speed up to compensate for the extra energy requirement due 
to a decrease in coupling energy, as shown in Fig. ^|(a), in which a large carrier velocity was chosen to exaggerate 
the effect. However, the behavior is richer than the one generally found for topological solitons because, in addition 
to the time dependent position Q(t), the ansatz contains an internal degree of freedom, £, so that the energy can 
also be transferred between different collective coordinates. With n < the last term of Eq. (|2^) decreases in the 
region of the defect, therefore f 2 has to increase accordingly. If n > and the breather is initially inside the defect, 
it simply slips away as in Fig. |2j(b) . If it is initially outside, for some suitable range of amplitude it is first slowed 
down and eventually reflected, as shown in Fig. ^(c); reflection occurs beyond r\ ~ 0.25. For smaller 77, breathers 
pass through the defect, indicating that the broader breathers are less influenced by the presence of defects, just as a 
large-wheel bike will not be stopped by a pebble or a ditch. In Fig.||(c) the breather has actually penetrated into the 
defect before being reflected. When the breather is trapped it oscillates between two positions, which may not be the 
defect boundary, as shown in Fig. 0(d). For values of 77 close to the threshold between trapping and non-trapping, 
the breather slowly turns around at the boundaries as in Fig. ||(e) . 

A necessary condition for a moving breather to be trapped in the above defect is k < 0. This statement can be 
proved through the following argument: a necessary condition for trapping is Q = more than twice, which, according 
to Eq. is equivalent to 

cosh 2 ( = 1 — cosh 2 (?7^) — Ksinh(rjl) cosh(rjl). (40) 

Since cosh 2 > 1, ko > and rj > 0, k has to be less than zero. We have therefore showed that trapping occurs only if 
the perturbed coupling constant is less than the unperturbed one, which is consistent with our simulations although 
it has been proven only in the collective coordinate approach. 

Since Eq. ( plO| ) contains only £, rj, a, and I, if the characters of the defect, i.e. the length, I, and the strength, k, have 
been fixed for a given system, and the initial position of the breather is chosen, the only factor which characterizes 
trapping is the breather amplitude. The initial value of <f> seems to have no consequence on the results. In general, if 
one finds that a breather passes through a defect for k < as in Fig.||(a), one can obtain trapping by increasing its 
amplitude. 

Because of the helicoidal structure of DNA, a given strand is alternatively inside and outside the bend so that it 
experiences a periodical modulation of its elasticity by an attached enzyme. We examined the consequence of such a 
modification by considering the following coupling constant modulation: 

k = k[Q(x + I) - 2G(x) + 9(a; - 0] (41) 

for which the coupling is first increased by n and then decreased by the same amount if k > 0. It can be viewed as a 
step approximation of one period of a sinusoidal modulation. In general one finds nothing essentially new with this 
perturbation because it is only a superposition of two step defects. However, this perturbation is asymmetric in space. 
By changing the sign of k we can reverse the orientation of the perturbation with respect to an incoming breather. 
If k < 0, a breather starting from the left side of the defect encounters first the region where the coupling constant is 
decreased. Table I summarizes the behavior of breathers with various amplitudes and the two possible signs of k and 
£. For negative k and positive £, the breather is reflected for intermediate rj while for large enough rj it is trapped. If 
one switches to positive k there is still a range of rj values that produce reflection, but for large rj the breather passes 
through. If the breather starts from the side where coupling constant is decreased the trapping can still exist even if 
the initial position of the breather is far away from the defect, but the pass-through region disappears as expected. 
These results show that it is the first encounter which determines the trapping. However, for this case of a composite 
defect the collective coordinate calculation can, in some cases, lead to qualitatively wrong results. The full numerical 
calculation shown in FigJJ(a) indicate that the breather can be trapped even if it were coming from the higher side of 
the defect. This points out the limit of the collective coordinate method for successive perturbations of the breather. 
The first interaction of the breather with a perturbation appears to be qualitatively well described. But then the 
perturbed breather is not accurately described by the ansatz. Thus when it encounters a second perturbation (here 
the second step in coupling constant), the collective coordinate description fails to describe the interaction. 
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V. DIRECT NUMERICAL SIMULATIONS 



Since the last example has shown that the collective coordinates cannot provide a full description of the breather 
dynamics, it is necessary to check them against full numerical simulations of Eqs. (jsj) . Using the breather solution 
given by Eq. fl22] ) as an initial condition, and periodic boundary conditions, we integrate Eqs. (||) with a 4 th order 
Runge-Kutta scheme and a time step chosen to provide a conservation of energy to an accuracy better than 1CU 6 
over a full simulation. The calculations have been tested on different system sizes to make sure that the results are 
not modified by boundary effects. The ansatz ( |22| ) is not an exact solution of the full set of equations because the 
transformation to the NLS form involved several approximations, however, except for very discrete cases or large 
amplitude breathers, it provides a rather good solution far away from the defect. As long as the breather is far away 
from the defect, one generally notices only a small decay of the initial energy peak due to radiation. 

Full numerical calculations have the advantage of allowing radiation and breaking of a breather. Furthermore, 
although the collective coordinate method starts from the perturbed NLS which requires small u e and u c and hence 
small amplitude, the full numerical calculations do not have this restriction. In what follows we show both energy 
distribution and the breather amplitude. The energy distribution is more relevant to the opening of DNA chain while 
the breather amplitude allows a comparison with the results of the collective coordinate calculations. 

Fi g. p| (a) is a typical case for trapping at an equivalent amplitude rj = 0.19. The correspondence is made from 
Eq. fl24[) and Eq. (p6|). The threshold for trapping predicted by collective coordinate is higher (77 = 1.01 for a breather 
initially at x = 12 when both systems have the same dimensionless coupling constant). As noticed earlier, it is 
not surprising to find such a discrepancy because we have used a sharp perturbation that violates the condition 
(|io|). However the full simulations confirm the qualitative predictions of the collective coordinate calculations: small 
amplitude breathers are transmitted, while larger ones are trapped. Fig. |^(b) shows that the energy distribution 
around the breather gets sharper in the region of the defect. Therefore a negative perturbation, which tends to trap 
breathers, is also favorable for base-pair opening since it concentrates the energy of the incoming breathers in a narrow 
domain. This sharpening of the breather shape occurs when the breather is inside the perturbation domain, whether 
it will stay trapped or not. One can also find on the contrary that, if a breather meets a positive perturbation, its 
energy distribution broadens. This behavior is similar to that of a vortex in shallow water: the vortex becomes wider 
when it is in shallower water and thinner in deeper water. ]l6| In the amplitude plot of Fig. ^(b), when the breather 
reaches the boundary of the defect, one can see two small reflected waves. They were not included in the collective 
coordinate analysis, and their presence explains part of the quantitative discrepancy between the analytical approach 
and the full simulations. Sometimes one can also notice that the breather changes its oscillation frequency after the 
collision with the defect. 

The results of the full numerical simulations show that, although the collective coordinate analysis is able to predict 
qualitatively the main features, in particular the existence of a threshold for trapping when the breather amplitudes 
increases, it is quantitatively wrong. The same conclusion had been found for an isolated impurity There are 
several reasons for that. Firstly we do not know an appropriate ansatz for the original equations of motion ([!]) and 
we start from a perturbed NLS Lagrangian which is already approximate. Then we use an ansatz which is localized 
in space and does not allow for the breaking of the breather or the emission of reflected waves. And finally the 
calculation assumes a smooth evolution of the coupling constant while we later use a sharp variation to make the 
analytical calculation possible. In spite of all their weaknesses, the collective coordinate calculations are however 
useful to get an insight of the behavior of the breather in the presence of the defect or even draw general conclusions 
on the kind of defects that can trap energy as explained above. 

Another point of interest is the trapping of several breathers in the defect region which could really enhance the 
energy density locally and cause local openings in DNA. We show in Fig. ^ examples for trapping for two kinds of 
coupling constant shapes: Fig. ||(a) is an example when the breather comes from the higher side of a two-step defect 
but is trapped. However, unless in favorable conditions we seldom find that two breathers can be trapped inside the 
same perturbation. When the second breather gets trapped it often kicks out the first breather that was trapped 
before, as shown in Fig. f|(b). In other cases we noticed that, when a first breather is trapped in the defect, a second 
breather that would have been trapped if it were alone, is on the contrary reflected. Therefore, if one studies only the 
positions of the breathers during their first interactions with the extended defect, it seems that the defect will never 
collect more than the energy of one breather. This is in fact not true, but the complete phenomena require a more 
detailed analysis. It is interesting to study the evolution of the energy in the region of the defect versus time. An 
example is shown on Fig. [| In this case the first breather that interacts with the defect has an amplitude 77 = 0.2 
which is above the trapping threshold and the second one has an amplitude 77 = 0.1 below the threshold. As expected 
the first breather is trapped and oscillates around the defect. The second one passes through the defect region that 
contains the first breather. However, if one looks at the energy density on the three-dimensional plot of Fig. ||(a), 
one can notice a significant increase of energy density after the interaction of the second breather with the defect. 
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The reason is that the second breather is only partly transmitted. A large part of its energy is given to the trapped 
breather, i.e. it stays in the defect region. The same phenomenon occurs again when the second breather collides a 
second time with the trapped breather. Due to this complex process, the time evolution of the energy inside the defect 
region (Fig. ^|-c) is a complicated curve, but it is important to notice that it tends to grow, and never falls again to 
a small value, indicating that the multiple collision process does cause a concentration of energy in the defect region. 
The origin of this localization of energy does not lie in breather trapping but in breather interactions in the presence 
of a perturbation, and therefore it is not included the collective coordinate description of Sect. IV. The result is very 
reminiscent of a mechanism described recently for energy localization due to discreteness effects in nonlinear lattices 
. In both cases the collisions of breathers, perturbed either by a defect or by discreteness, cause energy transfers 
that, on average, favor the big excitation at the expense of the small one. We have checked that the mechanism is not 
restricted to a particular case. Fig. ^ shows another example in which 3 breathers with the same initial amplitude 
7] = 0.2 were sent to the defect. Although the details of the process are different, they lead to the same final result: 
breather interactions in the presence of the defect tend to favor the formation of a large amplitude breather that 
concentrates a large part of the energy of the three incoming breathers and is finally trapped at the defect site. Hence 
the energy in the defect region settles to a high value. Tests have been performed with various breather amplitudes, 
leading to the same general result. 

VI. CONCLUSION 

Using a simple DNA model we have modeled the effect of a transcription enzyme by an extended modification 
of the coupling constant along the strands. The results show that such a perturbation is more efficient than an 
isolated impurity to trap breathers, in particular because trapping can occur provided that the amplitude of the 
incoming breather exceeds a threshold instead of requiring breathers with a well defined frequency. This conclusion 
can be derived from collective coordinate calculations as well as from numerical integration of the full set of equations 
of motion, although the collective coordinate method overestimates the trapping threshold. One cannot expect 
quantitative results from the collective coordinate analysis because we have violated at least one basic assumption, 
Eq. (|l(]|), to allow the analytical calculations, but it gives insight into the physics, and in particular a necessary 
condition for breather trapping which is confirmed by the full simulations. 

We have also showed that energy exchanges between a first breather, already trapped, and other incoming breathers 
can lead to a concentration of energy in the region of the defect. 

One may wonder whether the results obtained above for specific perturbations are extendable to more realistic cases. 
Although it is difficult to give general answers to this question, one can get insights through numerical simulations 
of the full system. In real DNA one has D = 0.03eV, a = 4.45A" 1 , k x = 0.08eV/A 2 , m = 300 a.m.u. for AT base 
pairs, while for GC pairs we have D = 0.035eV, k\ = 0.104. |l0| This is equivalent to k n « 0.13 to 0.15. In this 
range of coupling large amplitude breathers are trapped by discreteness |lq] . The collective coordinate calculations 
suggest that the low amplitude breathers, which can move, will not be trapped by the 20 base-pair defect. Simulations 
show that it is not necessarily so. For instance a 20 base-pair defect with K n = 0.12 in a chain with K n = 0.15 can 
trap breathers of various amplitudes. The energy exchange mechanism in the presence of the defect, discussed above, 
interferes with discreteness effects that can have similar effects to localize energy |17| |. Therefore, although we have 
exhibited a mechanism which is active in a wider frequency range than an isolated defect, the calculations performed 
on a simple model are not be sufficient to draw a conclusion about its validity to describe the effect of an enzyme on 
DNA transcription. It may however deserve attention because of its greater efficiency compared to the case of a point 
defect that was considered previously. 

In this work we have modeled the role of the enzyme by modulating only the coupling constant along the strands. 
As mentioned in the introduction, other possibilities could be considered, particularly if one attempts to take into 
account the enzyme specificity which suggests that the enzyme could have another role than merely bending locally 
the molecule. As a first step in this direction, we have considered a local change of the Morse potential in addition 
to the effect of the bending. Figure [7] shows the result of a numerical simulation where all the conditions are the 
same as for Fig. ||, except that, in addition to changing the coupling constant inside the defect to model the bending, 
we have also multiplied the denaturation energy of the base pairs (parameter D of Eq. (Q) ) by a factor 0.8. This 
means that we also assume that the enzyme can have some chemical effect to reduce the base pairing interaction. 
The comparison of Fig. 7 and Fig. 6c shows that this modification has a rather drastic effect on the results. It is 
easy to understand qualitatively why because locally the vibrating frequency of base pairs has been reduced. As we 
consider low energy breathers which are the most likely to be excited at physiological temperatures, their frequency, 
situated below the base-pair linear frequency of the unperturbed region because of the soft nonlinearity of the Morse 
potential, is however very close to the bottom of the phonon band of the unperturbed part of the molecule. As the 
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enzyme lowers the frequencies of phonon band in the defect region, the breather frequency is now in resonance with 
some modes of the phonon band of the defect. Therefore when the breather is trapped at the defect site by the 
bending, it is trapped in a region where it resonates with phonons. As a result it loses energy by radiation, but, as 
the emitted modes have a frequency below the lowest frequency of the unperturbed lattice, the vibrations are trapped 
in the defect region. One observes that the trapped breather spreads out its energy inside the defect region. When a 
second breather comes to this excited defect it is no longer repelled by a highly localized breather as in Fig. ^. Thus 
it is more likely to penetrate in the defect region too. This makes the energy localization effect more efficient and 
instead of the large oscillations that were observed in Fig. ^-c, Fig. shows that the energy in the defect region now 
grows steadily, each new breather having a high probability to add its contribution. Although it is still preliminary, 
this example shows that, if one combines the bending effect of the enzyme with some model for its specific action on 
the promoter site, one can perhaps provide a mechanism to achieve the local opening of the double helix which is 
required by DNA transcription. 
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FIG. 1. A one-dimensional lattice model of DNA. 
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FIG. 2. Evolution of collective coordinates in the presence of a step defect, Eq. ( poj ) for I = 10, (a) r\ — 0.5, k = —0.4, 
0(0) = 0.05, ^(O) = -0.06, C(0) = 15. (b) n = 0.21, k = 0.1, 0(0) = 0.4, 5(0) = 0.2, C(0) = -8. (c) n = 0.25, k = 0.5, 
0(0) = 0.4, ^(O) = -0.1, C(0) = 20. (d) 77 = 2.3, k = -0.5, 0(0) = 0.1, 5(0) = -0.13, C(0) = 25. (e) r) = 2.253, k = -0.5, 
0(0) = 0.1, 5(0) = -0.13, C(0) = 25. 
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FIG. 3. Half-tone plots of the energy distribution of breather evolutions, left, and the corresponding amplitudes, right, from 
direct numerical integration of Eq. (g) . In each figure the top insert shows the variation of the coupling constant used for the 
calculation, while the right inserts show snap-shots of the breather energy or amplitude distributions. Defects positions are 
shown in the plots axis, (a) K n = 0.44 outside the defect, K n = 0.22 for the defect. X(0) = 210, u c = 0.13, u e — —0.1. (b) 
K n = 1.04 outside the defect, K n = 0.52 for the defect. X(0) = 210, u c = 0.13, u e = -0.1. 
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FIG. 4. Half-tone plots of the energy distribution of two-breather evolutions, left, and the corresponding amplitudes, right, 
from direct numerical integration of Eq. (5) and various shapes of defects. In each figure the top insert shows the variation of 
coupling constant used for the calculation, while the inserts show snap-shots of the breather energy or amplitude distributions. 
The breathers start from X(0) = 120,230, defects positions are shown in the plots axis and (a) K n = 0.4 outside the defect, 
K n = 0.2 inside the defect and u c = 0.13, 0.11, u e = -0.11, -0.07. (b) K„ = 0.4 outside the defect, K n = 0.2 sin(n7r/Z) for the 
defect with n the distance from the center of the defect and I the half-length of the defect, and u c = 0.13, 0.12, u e — —0.11, —0.1. 
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FIG. 5. Interaction of two breathers with an extended defect. K„ = 1.0 outside the defect, K„ — 0.5 inside the de- 
fect which is 20 cells wide (between cells 145 and 165 in a 200-cell lattice) The parameters of the breathers are r\ = 0.1 
(tie = —0.014, u c = 0.333) and r\ = 0.2 (u e = —0.014, u c — 1.356), and both have a carrier wavector q = 0.1. (a) 
Three-dimensional picture of the energy density e n versus time, and (b) contour plot of the same energy density, (c) En- 
ergy E m inside the defect window (full line) and total energy of the chain E (dotted line) versus time. 
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FIG. 6. Interaction of three breathers with an amplitude r\ — 0.2 with the same extended defect as in Fig. 5. (a) Contour 
plot of the energy density in the chain versus time, (b) Energy E w inside the defect window (full line) and total energy of the 
chain E (dotted line) versus time. 
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FIG. 7. Energy in the defect region versus time for the interaction of three breathers with an extended defect. The conditions 
are exactly the same as in Fig. 6, but, in addition to changing locally the coupling constant inside the defect, we also reduce 
the barrier of the Morse potential by a factor 0.8 in this region to model some enzyme specific action. 
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TABLE I. Behavior of a breather determined by collective coordinate method with initial values (f> = —0.4, C, = — 13 and a 
two-step defect with length I = 10. The table lists the outcome of the interaction of the breather with the defect for different 
values of rj, and two different signs for k. 
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